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Abstract 

We consider a constrained minimax single facility location prob- 
lem on the plane with rectilinear distance. The feasible set of location 
points is restricted to rectangles with sides oriented at a 45 degrees 
angle to the axes of Cartesian coordinates. To solve the problem, an 
algebraic approach based on an extremal property of eigenvalues of 
irreducible matrices in idempotent algebra is applied. A new algebraic 
solution is given that reduces the problem to finding eigenvalues and 
eigenvectors of appropriately defined matrices. 

Key- Words: minimax location problem, rectilinear metric, idem- 
potent semificld, eigenvalues and eigenvectors. 



1 Introduction 

Models and methods of idempotent algebra [H EJ [31 [U [5] find increasing 
application in solving actual problems in engineering, economics, manage- 
ment, and other fields. Expressed in terms of idempotent algebra, a range of 
problems that are nonlinear in the ordinary sense, become linear, resulting 
in more simple analysis and solution techniques. Specifically, many classical 
optimization problems (in particular, in graph optimization and dynamic 
programming) reduce to solving linear vector equations, finding eigenval- 
ues and eigenvectors of matrices, and to working out similar computational 
procedures. 

Among the problems that have solutions in terms of idempotent alge- 
bra are certain location problems. A single facility one-dimensional location 
problem on a graph is examined in [61 [I] , where it is turned into a problem of 
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minimizing a rational function in the idempotent algebra sense. Note, how- 
ever, that the proposed solution concentrates on the analysis of polynomial 
and rational functions of one variable, and becomes less or no applicable in 
the case of two or more variables. 

In 0[8], a multidimensional constrained location problem on a graph is 
reduced to an optimization problem with a max-separable objective function 
that can be represented as a maximum of functions each depending only 
on one variable. An efficient computational procedure is proposed which, 
however, seems to have limited application only to location problems with 
max-separable objective functions. 

In this paper, we consider a constrained rectilinear minimax single facil- 
ity location problem on the plane, where the feasible set of location points 
is restricted to rectangles that have sides oriented at 45° angle to the axes 
of Cartesian coordinates. In the unconstrained case, the problem has a well- 
known closed- form solution derived in [£l [10] from geometrical arguments. 

A new algebraic solution to the unconstrained problem is given in [11] 
on the basis of the approach proposed in |12t [TT] to solve minimax location 
problems with rectilinear and Chebyshev metrics. The approach exploits 
properties of eigenvalues of irreducible matrices in idempotent algebra to 
reduce certain location problems to finding eigenvalues and eigenvectors of 
appropriately defined matrices. The aim of the paper is to extend the alge- 
braic solution to the constrained location problem under consideration. 

We start with a brief overview of basic facts about idempotent algebra 
to provide a background for subsequent sections. Some relevant results from 
the spectral theory of matrices are outlined including an extremal property 
of eigenvalues of irreducible matrices. We provide a conventional descrip- 
tion for the location problem of interest and then represent the problem 
in terms of idempotent algebra. An algebraic solution is first given to the 
unconstrained problem and then extended to the constrained case. 

2 Preliminary Results 

In this section, a short introduction to idempotent algebra based on [13\ [T3] 
is given as background material to the solution approach. Further detail on 
the idempotent algebra theory and methods is provided in [IJ [2j O HI [5] . 

2.1 Idempotent Semifield 

Consider a set X endowed with two operations, addition © and multipli- 
cation ®. We assume (X, 0, 1, ©, (g>) to be a commutative semiring with 
zero and identity 1, where addition is idempotent and multiplication is 
invertible. Since the nonzero elements in X form a group with respect to 
multiplication, this semiring is often referred to as the idempotent semifield. 
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The power notation is denned as usual. Let X + = X \ {0}. For each 
x € X + and any integer p > 0, we have x — 1, (D p = 0, x p = x^ 1 <g> x = 
x <8) x p ~ x , and x~ p = (x" 1 ) 10 . 

We assume that the integer power can naturally be extended to the cases 
of rational and real exponents. 

In the algebraic expressions below, the multiplication sign © , as usual, is 
omitted. The power notation is thought in the sense of idempotent algebra. 

The idempotent addition induces a partial order on X so that x < y if 
and only if x@y = y. With this definition, it is easy to verify that x < x(By, 
y < x © y , and that both addition and multiplication are isotonic. 

Finally, it is not difficult to show that the binomial identity takes the 
form (x © y) a = x a © y a for all a > 0. 

As an example of such semirings, one can consider the idempotent semi- 
field of real numbers 

Kmax,+ = {RU {-oo}, -oo, 0, max, +). 

In R max ,+ , the null element is — oo, and the identity is 0. For each 
x € K, there exists the inverse x _1 that is equal to — x in the conventional 
arithmetic. For any x,y € R, the notation x y conforms to the arithmetic 
product xy . The partial order coincides with the natural linear order in R . 

2.2 Vector and Matrix Algebra 

Vector and matrix operations are routinely introduced on the basis of the 
scalar operations defined on X. Consider a Cartesian product X™ with 
its elements represented as column vectors. For any vectors x = (xj) and 
y = (yi), and a scalar c, vector addition and multiplication by scalars follow 
the conventional rules 

{x © y}i = Xj © yi, {cx}i = cxi. 

As usual, a vector y £ X n is linearly dependent on vectors x±, . . . , x m £ 
X™ if there are scalars ci, . . . , c m G X such that y = c\X\ © • • • © c m x m . In 
particular, y is collinear with x when y = ex. 

For any column vector x = (x«) € X™ , we define a row vector x = (x~) 
with elements x^ = x~ l . For all x, y € X™ , the component-wise inequality 
x < y leads to the inequality x~ > y . 

For any conforming matrices A = (a^), B = (bij), and C = (cij) with 
entries in X , matrix addition and multiplication together with multiplication 
by a scalar c are performed in accordance with the formulas 

{^4 © B}ij = dij © bij, {BC}ij = bikCkj, {cA}ij = cdij. 

k 

A matrix with all entries equal to zero is called the zero matrix and 
denoted by 0. 
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As is customary, a square matrix that has the diagonal entries equal to 
1, and all off-diagonal entries equal to is called the identity matrix and 
denoted by /. 

For any square matrix A ^ and integer p > 0, we have A = I, 
A p = ap^A = AA p ~ l . The trace of a matrix A = ( aij ) € X nxn is given by 
tr A = an © • • • © a nn . 

A matrix is irreducible if it cannot be put into a block triangular form 
by simultaneous permutations of its rows and collumns. 

2.3 Eigenvalues and Eigenvectors 

A scalar A is an eigenvalue of a matrix A 6 x nxn if there exists a vector 
x € X™ \ {0} to provide the equality Ax = Xx. Any vector x that satisfies 
this equality is an eigenvector of A with the eigenvalue A. 
Any irreducible matrix A has one eigenvalue given by 

n 

X = Q)tr 1 / k (A k ). (1) 

k=l 

The eigenvectors of the matrix A have no zero entries. 
For any irreducible matrix A £ X nxn , the eigenvectors corresponding to 
A can be found as follows. First we evaluate the matrix 

A x = X~ X A © • • • © X~ n A n . 

Let af be column i in A x , and o 4 x its diagonal element. We take the 
set of all columns af such that =1, and then reduce it by removing 
those elements, if any, that are linearly dependent on others. Finally, the 
columns in the reduced set are put together to form a matrix A + . 

The set of eigenvectors that corresponds to A (including the zero vector) 
coincides with the linear span of the columns of A + , whereas any eigenvector 
is represented as 

x = A + s, 

where s is a nonzero vector of appropriate size. 

We conclude with the results from [Til 02] that underlies the solutions 
in the subsequent sections. 

Lemma 1. Let A = (aij) G x nxn be an irreducible matrix with an eigen- 
value A. Suppose u = (ui) and v = (vi) are eigenvectors of the matrices A 
and A T . 

Then it holds that 

min x Ax = A, 
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where the minimum is attained at any vector 

( « _1 



x 



< a < 1. 



(2) 



a„,a— 1 



3 A Minimax Location Problem 

In this section, we give two equivalent representations of the constrained 
minimax single facility location problem on a plane with rectilinear distance. 
We start with a problem definition in the usual notation, and then suggest 
an alternative form written in terms of the semifield R m ax,+ • 

3.1 Conventional Representation 

The distance between any two vectors r = (ri,r2) T and s = (si,s 2 ) T on 
the plane R 2 in rectilinear metric is 

p(r,s) = |n - si| + |r 2 - s 2 |. (3) 

Suppose there are vectors T{ = (Vij,r 2 j) T € R 2 and numbers Wi £ R 
for i = 1, . . . , m, where m > 2. Given a feasible set S C R 2 , the location 
problem of interest is to find the vectors x = (xi,X2) T that provide for 

min max (/o(r*j, x) + Wi). (4) 

xdS l<i<m 

In this paper, we restrict our consideration to feasible location sets that 
have the form of rectangles with sides oriented at 45° angle to the axes of 
Cartesian coordinates. 



3.2 Algebraic Representation 

First we represent the metric ([3|) in terms of operations in the semifield 

Rmax,+ 

p(r, s) = (s^Vi © r 1 _1 si)(s 2 _1 r2 © r^ 1 s 2 ). 
Now problem can be rewritten in the form 

mini^(x), (5) 
xeS 

where the objective function is given by 

m 

¥( x ) = ® 1 r u © rf- 1 xi)(x^ 1 r 2i © r 2i 1 x 2 ). 

i=l 

To solve this constrained location problem we first examine its uncon- 
strained version. 
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4 The Unconstrained Problem 



Let us consider problem Q and put 5 = R 2 to get the problem 

min <p(x), (6) 
xeR 2 

For this unconstrained problem, a solution based on geometrical argu- 
ments is obtained in [9] , where the problem is called the rectilinear messenger 
boy problem. A similar solution in the case when Wi = for all i is given 
in [10]. 

Now we outline how the problem can be rearranged and then solved 
algebraically. Additional detail is given in 



4.1 Further Algebraic Transformation 

Consider the objective function of the problem. After some simple algebra, 
we get 

<p(x) = ax± 1 X2 © bx±X2 1 © cx± 1 X2 1 © dxiX2, (7) 



where 



a = Qw i r li r 2i 1 , b = (J) w i r li 1 r 2i , (8) 
i=i i=i 

m m 

0^ur2i, rf = u; i r ri lr 2i 1 - ( 9 ) 



c = 

8=1 8=1 



Let us introduce a vector and a matrix A of third order as follows 

/ x x \ / a \ 

y = x 2 , A = 6 c . (10) 
\ x^ 1 / \ d / 

Now the objective function can be represented as 

ip(x) = y Ay, 

whereas problem ([!]) reduces to that of the form 

min y~ Ay. (11) 

It is easy to see that the solutions of the new extended three-dimensional 
problem do not all correspond to solutions of ©. By virtue of the definitions 
at (fTOj) , a proper solution of the extended problem is a vector y such that 
its first and last components are inverse to each other. 
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4.2 Algebraic Solution 



Let us examine extended problem (jlip . First note that the matrix A defined 
in (|10p is irreducible. From Lemma [1] it follows that the minimum in (jllf) is 
equal to the eigenvalue A of A , and the minimum is attained at any vector 
defined according to ([2]) through eigenvectors of the matrices A and A T . 
Application of ([U gives 

A = (a6 0cd) 1/2 . 

Now we obtain the corresponding eigenvectors of the matrices A and 
A T . Note that A T is derived from A by interchanging a and b, as well as 
c and d. Therefore, we can first find the eigenvector of A, and then turn it 
into the eigenvector of A T by the above interchange. 

To get the eigenvectors of A, we need to examine the matrix 

3 / X~ 2 ab A _1 a A~ 2 ac \ 

A x =Q)\- k A k = \- l b 1 \~ l c . 
fc=l V \- 2 bd X~ x d \~ 2 cd ) 

There are three cases to investigate. 



4.2.1 The case A 2 = ab > cd 

The matrix A x takes the form 

/ 1 A _1 a A~ 2 ac 
A x = X^b 1 A^c 
\ \-%d X^d X~ 2 cd 

Since X~ 2 cd < 1 , only the first two columns of A x are eigenvectors of A. 
Considering that both columns are collinear, we take only one of them, say 
the first column, to form the matrix A + . Now we arrive at the eigenvectors 
of A defined as 

1 

u = A + s = | A^fr |s, s G R. 

A~ 2 M 

By interchanging a and b, and then c and d, we get the eigenvectors of 

1 

X~ l a | t, t€R. 
X~ 2 ac 

Application of ([2]) leads to the solution of the extended problem in the 
form 

1 

y | A 1 - 2 "^- 1 ^ | s a t a ~\ s,ieR. 

X 2(l-2a) aa -l b a c a~l d a 
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Under the condition that the first and last components of the vector y 
must be reciprocal, we have 

s^' 1 = (a a b a - l c l - a d- a ) 1/2 . 

Turning back to problem we arrive at 



x 



(a^-V^cT ) 1 / 2 
(a a - 1 b a c 1 - a d-° 1 ) 1 / 2 



4.2.2 The case A 2 = cd > ab 
In a similar way as above, we get the matrices 

A x = 

The eigenvectors of A and A T take the form 



A~ 


2 ab 


X~ l a 


X~ 2 ac \ 






A 


^b 


1 


X~ l c 






A" 


2 bd 


\~ l d 


1 J 




V i / 



u 



( ^ 




( x- 2 bd \ 






\s, v=\ 


X^d 


\t, s,t€ 


\ t J 




\ 1 J 





The extended problem has the solution 

/ A 2 ( 1 - 2a )a a 6 a - 1 c Q d a - 1 \ 
y=i A 1 - 2 "^^" 1 I s a t a -\ s,teR. 

To get a solution of ([6|), we put 

s^- 1 = {a-^^c^d 01 ) 1 / 2 , 

and then get 



x 



(a- a b l - a c a d a - 1 ) 1 / 2 



4.2.3 The case A 2 = ab = cd 

From the matrix A x we successively get 

1 



i / (a a 6 Q - 1 c 1 - Q d- a ) 1 / 2 

A o , a; - 1 
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4.3 Summary of Solutions 

It is not difficult to verify that the above three solutions can be combined 
into one to represent the result as follows. 

Lemma 2. With the notation ([8]) the minimum in ([6]) is equal to 

A = (abecd) 1 / 2 

and attained at the vector 



x 



(a a b a - 1 c 1 ~ a d- a ) 1 / 2 

X 2a-l( a -a b l-a c l-a d -ay/2 



for all a such that < a < 1 . 

In terms of ordinary arithmetic operations, we have the following state- 
ment. 

Corollary 3. Suppose that 

a= max (wi + r u - r 2 i), b= max (wi - m + r 2 j), 

l<i<m l<i<m 

c = max (uui + m + r 2 i), d= max (w{ - m - r 2 j). 

l<i<m l<i<m 

Then the minimum in is given by 

A = max(a + b, c + d)/2, 
and it is attained at the vector 



x 



f (o-d) - ^(b-c) 



(2a - 1)A - f (a + d) + ±=^(b + c) 
for all a such that < a < 1 . 

5 Constrained Problems 

The main idea that underlies the solution of the constrained problem under 
consideration is to put it in an unconstrained form, and then take advantage 
of results in the previous section. 

First we take the unconstrained problem at ©, and represent it in a 
slightly different way that is more convenient for further incorporation of 
constraints. With subscripts added to the above notation, we have 

m m 

°o = ® Wirur^j 1 , b = (J) Wir^r 2i , 

i=l i=l 

m m 

c o = (J) w i r ur2i, d = (B^ifT/r" 1 , 

i=l i=l 
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and 



A = (a 6 © c do) 7 • 
Furthermore, we define an objective function 

<Pq(x) = Aq 1 (a x^ 1 x 2 © boxix^ 1 © cqx^x^ 1 © d xix 2 ). 



Now consider the problem 



xeR 2 



min LfQ (x). 



It is easy to see that this problem is equivalent to ([6]) in the sense that 
they both have the same solution set. Note, however, that the minimum in 
the last problem is 1 = 0, whereas it is equal to A = Ao in ([6]). 

5.1 Location in 45° Rotated Rectangles 

Given numbers a%, b\, ci, d\ , we assume the feasible location set S to be the 
intersection of the half-planes defined in standard notation by the inequali- 



It is easy to see that the intersection of the half-planes, if nonempty, has 
the form of an upright rectangle rotated by 45° around its center. 

In terms of the semifield K m ax,+ , the inequalities can be rewritten as 



ties 



x\ + X2 + a\ < 

■X\ — X2 + C\ < 



x\ - x 2 + h < 0, 
x\ + x 2 + d\ < 0. 




b\Xix 2 < 1 
dixix 2 < 1 



and then further combined into one inequality 



a\x 1 1 x 2 © b\Xix 2 1 © c\x x l x 2 1 © d\X\x 2 < t. 



Let us introduce a function 



ifi(x) = aix 1 x 2 © bix±x 2 © c\X x x 2 © d\X\x 2 

and note that 9?i(a^) < 1 if and only if x G S . 
We define an objective function 



V>(x) = ipo{x) © ipi{x). 



Furthermore, with the notation 



a = A ao © ai 
c = \q 1 cq © ci, 



b = \ %®bi, 
d = Aq 1 ^ © d\, 
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we can put the objective function in the form of ([7]) 

ifj(x) = ax 1 1 x 2 © bx\%2 1 © cx^ x 2 1 © dX\X 2 . 
Now the constrained problem ([5]) reduces to an unconstrained problem 

min ip(x) (12) 

x£R 2 

that has the form of Q and the solution given by Lemma 
5.2 Solution to the Constrained Problem 

First note that the minimum in problem (|12j) is given by A = (a&© cd) 1 ^ 2 > 
1 . It is not difficult to see that A = 1 if and only if the feasible set of the 
original constrained problem is nonempty and it has an intersection with 
the solution of the associated unconstrained problem. 

In the case of A > 1, the solution to (fT2"j) can be considered as an 
approximate solution for the constrained problem, that is of independent 
interest. 

In the conventional notation, the solution to (|12p can be described as 
follows. 

Corollary 4. Suppose that 

a = max(ao — Ao, «i), b = max(&o — Ao, 

c = max(co — Ao, ci), d = max(do — Ao, di), 

where 

a = max (im + r u - r 2 i), 

l<i<m 

b = max (wi - r u + r 2 {), 

l<i<m 

c = max (im + ru + r 2 i), 

l<i<m 

d = max (im - ru - r 2 i), 

l<i<m 

A = max(a + ^o, c + <i )/2. 
Then the minimum in (|12p is given by 

A = max(a + b, c + d)/2, 
and it is attained at the vector 




f(a-d) - ^(b-c) 
(2a - 1)A - f (a + d) + + c) 



for all a such that < a < 1 . 
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